Everything is Entangled 
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We show that big bang cosmology imphes a high degree of entanglement of particles in the 
universe. In fact, a typical particle is entangled with many particles far outside our horizon. However, 
the entanglement is spread nearly uniformly so that two randomly chosen particles are unlikely to 
be directly entangled with each other - the reduced density matrix describing any pair is likely to 
be separable. 



ERGODICITY AND PROPERTIES OF TYPICAL 
PURE STATES 



When two particles interact, their quantum states gen- 
erally become entangled. Further interaction with other 
particles spreads the entanglement far and wide. Subse- 
quent local manipulations of separated particles cannot, 
in the absence of quantum communication, undo the en- 
tanglement. We know from big bang cosmology that our 
universe was in thermal equilibrium at early times, and 
we believe, due to the uniformity of the cosmic microwave 
background, that regions which today are out of causal 
contact were once in equilibrium with each other. Be- 
low we show that these simple observations allow us to 
characterize many aspects of cosmological entanglement. 

We will utilize the properties of typical pure states in 
quantum mechanics. These are states which dominate 
the Hilbert measure. The ergodic theorem proved by 
von Neumann [l[ implies that under Schrodinger evolu- 
tion most systems spend almost all their time in typi- 
cal states. Indeed, systems in thermal equilibrium have 
nearly maximal entropy and hence must be typical. Typ- 
ical states are maximally entangled (see below) and the 
approach to equilibrium can be thought of in terms of 
the spread of entanglement. 

Consider a large system subject to a linear constraint 
R (for example, that it be in a superposition of energy 
eigenstates with the energy eigenvalues all being near 
some E^), which reduces its Hilbert space from TL to 
a subspace Hr. Divide the system into a subsystem A, 
to be measured, and the remaining degrees of freedom 
which constitute an environment B, so H = Ha (Xi Hb 
and 



PA = PaW = TrB|-0)(V'| 



(1) 



is the density matrix which governs measurements on A 
for a given pure state ip of the whole system. Note the as- 
sumption that these measurements are local to A, hence 
the trace over B. (Sec Fig. [TJ) 

It can be shown ^13|] (see also [1,[5[), using the concen- 
tration of measure on hyperspheres 6| (Levy's Lemma) , 




FIG. 1: The entire system is in a pure state t/; subject to a 
constraint on total energy. Tracing over the shaded region B 
yields a density matrix pA- For typical tp (which dominate 
the set of possible states), pA is nearly thermal and nearly 
maximally entangled. 



that for almost all tfj G Hr, 



Pa(V') ~ Trs (p*) = Qa 



(2) 



where = Ir/cIr is the equiprobable maximally mixed 
state on the restricted Hilbert space %r {1r is the iden- 
tity projection on T-Lr and d^, the dimensionality oH-Lr). 
U,A = Trs (p*) is the corresponding canonical state of 
the subsystem A. The result holds as long as ds ^ d^, 
where cIa and cLb are the dimensionalities of the I-La and 
T-Lb Hilbert spaces. (Recall that these dimensionalities 
grow exponentially with the number of degrees of free- 
dom. The Hilbert space of an n qubit system is 2"- 
dimensional.) In the case of an energy constraint R, U,a 
describes a perfectly thermalized subsystem with tem- 
perature determined by the total energy of the system 
(i.e., a micro canonical ensemble). 

To state the theorem in [2] more precisely, the 
(measurement-theoretic) notion of the trace-norm is re- 
quired, which can be used to characterize the distance 
between two mixed states pA and U,a'- 



\PA 



JIaIIi = Try' [pA - ^Af 



(3) 



This quantifies how easily the two states can be distin- 
guished by measurements, according to the identity 

IIPA-f^Alli -sup||o||<iTr(p^O-r!^0) , (4) 
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where the supremum runs over all observables O with op- 
erator norm ||0|| < 1. The trace on the right-hand side of 
(HI is the difference of the observable averages (O) evalu- 
ated on the two states pA and Q.a, and therefore specifies 
the experimental accuracy necessary to distinguish these 
states in measurements of O. 

The theorem then states that (for e > 0) 



value hvdA'- 



Prob 



Ipa (V') - Trs (p*) 111 > e + d^d^^^^ 
< 2exp(-e2dfl/187r3) . 



(5) 



In words: let ip be chosen randomly (according to the 
Haar measure on the Hilbert space) out of the space of 
allowed states T-Lr] the probability that a measurement 
on the subsystem A only, with measurement accuracy 
given by the rhs of (O, will be able to tell the pure state 
■0 (of the entire system) apart from the maximally mixed 
state p* is exponentially small in d/j, the dimension of 
the space Hr of allowed states. Conversely, for almost 
all pure states tp any small subsystem A will be found 
to be extremely close to perfectly thermalized (assum- 
ing the constraint R on the whole system was an energy 
constraint). 

As mentioned, the overwhelming dominance of typi- 
cal states ip is due to the geometry of high-dimensional 
Hilbert space and the resulting concentration of measure. 
It is a consequence of kinematics only - no assumptions 
have been made about the dynamics. Almost any dy- 
namics - i.e., choice of Hamiltonian and resulting uni- 
tary evolution of ip - leads to the system spending nearly 
all of its time in typical states for which the density ma- 
trix describing any small subsystem A is nearly thermal 
[II) 0| ■ Typical states ip are maximally entangled, and the 
approach to equilibrium can be thought of in terms of the 
spread of entanglement, as opposed to the more familiar 
non-equilibrium kinetic equations. 

Since generic pure states tend to evolve into typical 
states, any mixture of pure states is likely to evolve into 
a mixture of typical states. Hence, our analysis does 
not require any specific assumptions about whether the 
system (i.e., the universe) is in a pure or mixed state. If 
it is in a mixture, we simply have (classical) probabilities 
of finding the system in one of two or more typical pure 
states. For simplicity, in the rest of the paper we will 
always assume the system as a whole is in a pure state. 

We can restate these results in terms of the entangle- 
ment entropy of the subsystem A, thereby making con- 
tact with the Second Law of Thermodynamics. The en- 
tanglement entropy is simply the von Neumann entropy 
of pa: 



S{pa) = -Tr PA In PA 



(6) 



Using the same results on the concentration of measure, 
it can be shown [8j that, for the overwhelming majority 
of pure states ip, S{pa) is extremely close to its maximum 



a 



Prob [S'(pa) < IndA 

{dAde - l)Ca^ 



< exp 



(lndA)2 



(7) 



where a > 0, /3 = dA/ds and C = (Stt^)"!. This im- 
plies [l|, 0] that, for almost any choice of dynamics, a 
subsystem A is overwhelmingly likely to be found with 
nearly maximal entropy S{pa)- The Second Law is seen 
to hold, in a probabilistic sense, even though the underly- 
ing dynamics is time- reversal invariant: density matrices 
Pa with small entropy are highly improbable, and if A is 
found in a low-entropy state, the entropy is more likely 
to increase than decrease over any macroscopic time in- 
terval. 



COSMOLOGY 

In the following discussion we assume a semiclassical 
space-time and focus on field-theoretic degrees of freedom 
(e.g., particles such as photons or electrons). Although 
the analysis takes place in curved space, quantum grav- 
itational effects are never significant, and the rules of 
ordinary quantum mechanics apply throughout (just as 
they do in the Earth's gravitational field). We adopt a 
cosmological time coordinate (e.g., that of the FRW met- 
ric) and evolve the collective wave function of particles 
using the Schrodinger equation in those coordinates. 

The cosmic microwave background provides direct ev- 
idence that the universe was in thermal equilibrium at 
temperatures of order eV. Nucleosynthesis of light ele- 
ments provides indirect evidence of thermal equilibrium 
at temperatures in the keV to MeV range. This sug- 
gests that the state describing the universe in the past 
was typical. The ergodic theorem [l|, or equivalently, the 
Second Law, implies that the universe is likely to be in 
a typical state today. Thus the entanglement entropy of 
any subsystem A is likely close to maximal. 

To proceed further we recall that a cosmological hori- 
zon volume is the largest region over which causal contact 
is possible. The size of this region is 



dHit)^a{t) 



dt' 
^ 



(8) 



where a(t) is the FRW scale factor. Our present hori- 
zon volume consists of many sub-regions that are only 
now coming into causal contact, at least as implied by 
([5]) under ordinary (e.g., radiation- or matter-dominated) 
expansion. The fact that the temperature and distribu- 
tion of density perturbations (not to mention stars and 
galaxies) are approximately uniform over these regions 
suggests that, somehow, they were already in causal con- 
tact during some previous epoch. Most researchers now 
believe that this is due to a period of exponential growth 
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in a{t) known as inflation. During this era the metric 
was approximately that of de Sitter space and the en- 
ergy density was dominated by the vacuum energy of the 
inflaton field. In this scenario, the currently visible uni- 
verse originated from a progenitor region much smaller 
than the horizon volume at the start of inflation: 

^0 < rfff(iinflation) ■ (9) 

Because the entire horizon volume at t = ^inflation was in 
equilibrium, all of our visible universe and regions which 
will only become visible in billions of years experienced 
similar initial conditions, thus explaining the observed 
homogeneity and isotropy. At the quantum level, this 
equilibrium assumption implies that the pure state ip de- 
scribing a region of size d/f (^inflation) at the start of infla- 
tion was typical. Due to the inequality ([9|), entanglement 
today must extend far beyond the currently visible uni- 
verse. In fact, as we show below, particles in our horizon 
volume are mostly entangled with particles outside of it. 

We can see explicitly how entanglement is transferred 
by considering the inflaton field in the standard model of 
slow-roll inflation. Before inflation begins, the inflaton 
and other degrees of freedom are in thermal equilibrium 
and we expect their states to be typical. Once the infla- 
ton vacuum energy begins to dominate the stress energy 
tensor, the universe supercools and the gravitational dy- 
namics is determined by the semiclassical evolution of the 
scalar field as it slowly rolls along its nearly fiat potential. 
However, at the quantum level the inflaton wave function 
at each position in space is still entangled with the wave 
function at other positions: due to the non-zero de Sitter 
temperature and inflaton-graviton scattering, there are 
interactions which "measure" the value of the inflaton 
field and entangle its wave function with nearby degrees 
of freedom. (Indeed, this has to be the case for a semi- 
classical space-time to emerge, whose dynamics is mainly 
driven by the vacuum energy of the infiaton.) We need 
only require that the ergodic theorem apply to the sys- 
tem comprised of infiaton field, background particles and 
gravitons, during an epoch in which gravitational effects 
are small (e.g., the de Sitter timescale is larger than par- 
ticle interaction timescales). Under this assumption, the 
wave functions of the infiaton field in different patches of 
each horizon volume are entangled. 

Due to the exponential expansion in de Sitter space the 
physical separation between points grows superluminally 
and regions which were originally in the same horizon 
volume of size ^^(tinflation) become space- like separated. 
Local evolution in each causally separate region cannot 
undo the pre-existing entanglement. When the inflaton 
field decays, causing particle production and reheating 
the universe, this entanglement is transferred to the de- 
cay products, which include the particles that make up 
the universe today. The discussion above is in the con- 
text of a speciflc model of inflation, but in general we 
expect any mechanism for superluminal expansion which 



solves the isotropy and homogeneity problems will lead 
to entanglement across many horizon volumes. 

To summarize, modern cosmology suggests that most 
of the particles in the visible universe exhibit a high 
degree of entanglement with degrees of freedom far be- 
yond our horizon volume. While it is true that gravi- 
tational clumping (e.g., of galaxies or stars) Q allows 
local deviation from thermal equilibrium, entanglement 
with causally separated regions produced in earlier cos- 
mological epochs cannot be removed by subsequent local 
dynamics. 



ASPECTS OF MAXIMAL ENTANGLEMENT 

I. Maximal entanglement and Schmidt decomposition 
We argued above that any small subsystem A ( "small" 
here includes our entire horizon volume today!) is max- 
imally entangled with the rest of the universe (most of 
which is not yet visible to us). That is, S{A) fa \ndA- 
We can further interpret this using the Schmidt decom- 
position theorem: for any pure state iI^ab of a composite 
system AB, there exist orthonormal states ipA^^ for sys- 
tem A and V's'"' for system B such that 

^AB^Yl (10) 

n 

1/2 

where A„ are nonncgative real numbers satisfying 
X^n -^n = 1- This is a simple consequence of the sin- 
gular value decomposition theorem. The dimensions of 
"Ha and Hb can be very different, and the range over 
which the sum in Eq. (|10|) runs is determined by the 
smaller Hilbert space. Note that the Schmidt decompo- 
sition states might be quite complex - possibly involving 
superpositions of many degrees of freedom. Tracing over 
B yields a density matrix pA with eigenvalues A„. From 
the previous discussion we know that all A„ « 

A measurement of subsystem A which determines it 

(n) 

to be in state ip\ implies that the rest of the universe 

must be in state For example, A might consist of a 

few spins [13] ; it is interesting, and perhaps unexpected, 
that a measurement of these spins places the rest of the 
universe into a particular state V'b'''- As we will see be- 
low, in the cosmological context these modes are spread 
throughout the universe, mostly beyond our horizon. Be- 
cause we do not have access to these modes, they do not 
necessarily prevent us from detecting A in a superposi- 

(n) 

tion of two or more of the ip\ . 

However, if we had sufficient access to B degrees of 
freedom (for example, if the relevant information diffcr- 

(n) 

entiating between ipg states is readily accessible in our 
local environment or in our memory records) , then the A 
system would decohere into one of the . 



FIG. 2: Tracing over the shaded region B yields a density 
matrix pxY- If the dimensionality ds is much smaller than 
dxdy, the entanglement of formation Ef{pxY) will be nearly 
maximal. Taking B to be a horizon volume, this implies that 
X is mostly entangled with degrees of freedom outside its own 
horizon. 

II. Entanglement across horizons 

Consider Fig. ([2]), where BUX is our horizon volume, 
with X a small subregion. Tracing over B yields a den- 
sity matrix pxY describing the entanglement of region 
X with the rest of the universe Y (all of which is out- 
side the currently visible universe; Y is assumed much 
larger than B). Because entanglement should be roughly 
uniformly distributed over degrees of freedom in a typ- 
ical state, we expect that most of the entanglement en- 
tropy S{pxy) (which must be nearly maximal) is with 
modes in Y rather than B. Indeed, one can show (the- 
orem V.l of [1]) that if ds <C dxdy (i.e., many more 
degrees of freedom in X iJY than in _B), it is exponen- 
tially likely that the entanglement of formation Ef{XY) 
is close to log(c?xrfy) (i-e., is also nearly maximal). The 
entanglement of formation is a measure of entanglement 
for mixed states, such as pxv illl- It is defined as the 
minimum entanglement resource necessary to create pxY 
without further quantum communication. Alternatively, 
Ef{XY) is equal to the least expected entanglement of 
any ensemble of pure states which realize pxY- That is, 
for all decompositions 

Pxy=E^''l'^W)(0WI , (11) 

where Wi are probabilities and (/'xv ^ pure 
state, Ef{XY) is the minimum expected entanglement 
^ WiS'((/)3fy). These statements imply that most of the 
entanglement entropy S{pxy) is due to entanglement 
with modes of Y , which are causally disconnected (space- 
like separated) from X. 

III. Small systems are likely to he in separable states 
Fig. dS]) depicts two small regions X and Y (although 

depicted as far apart, they could also be spatially prox- 
imate). For example, each could consist of a single or 
small number of individual particles. The approximately 
uniform distribution of entanglement over all degrees of 
freedom in a typical state suggests that X and Y share 



FIG. 3: Tracing over the shaded region B yields a density 
matrix pxY- If X and Y are small subsystems {{dxdy)^ ^ 
ds), pxY is likely to be separable - X and Y are not directly 
entangled with each other. Note X and Y need not be widely 
separated; the figure is drawn this way for visual clarity. 



only negligible entanglement directly with each other. A 
measure of this direct entanglement is again the entangle- 
ment of formation for the density matrix pxY '■ Ef{XY), 
which we expect to be small. Indeed, theorem V.l in 
8*] provides an upper bound on Ef(XY) which vanishes 
in the limit that is much larger than dxdy- When 
ds S> (dxrfy)^, pxY is exponentially likely to be separa- 
ble: 

PXY ^J2^'P^'^Py ' (1^) 

where Wi are real, positive and sum to unity, and Px and 
Py are density matrices on X and Y. Separable states 
may exhibit classical correlations, but no entanglement. 
Even the classical correlations must be small because we 
know that pxY ~ ^xy /(dxdY)- 



CONCLUSIONS 

The cosmological quantum state is likely to be typical 
in a Hilbert space describing degrees of freedom over a 
region many times as large as the visible universe (our 
current horizon volume). This implies a high degree 
of entanglement of particles, with the entanglement dis- 
tributed uniformly over most of the degrees of freedom. 
As a consequence, small subsystems are mostly entangled 
with particles far beyond the horizon, and two randomly 
chosen small subsystems are unlikely to be directly en- 
tangled with each other. 
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